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RESUMEN 

Estudiamos la disipacion turbulenta en el MIE y algunas implicaciones para 
la formacion y evolucion de galaxias usando simulacioncs MHD numcricas de fluidos 
compresibles en 2D. La energia cinetica Ek se inyecta por fuentes estelares formadas 
autoconsistentemente. En el fluido coexisten el regimen de turbulencia forzada 
y el de decaimiento. En la regiones turbulentas activas, Ek es disipada local y 
ehcientemente. En el regimen de decaimiento E}.(t) decae ~ (1 + t)~°- s . Los 
movimientos turbulentos residuales pueden propagarse distancias del orden del alto 
del disco gaseoso lo cual sugiere que la turbulencia puede propiciar el soporte vertical 
con una formacion estelar autoregulada al nivel del disco, mas no al nivel de todo 
el halo cosmologico como se requiere en ciertos modelos de formacion de galaxias. 

ABSTRACT 

We study turbulent dissipation in the ISM and explore some implications for 
galaxy formation and evolution using 2D MHD numerical simulations of compress- 
ible fluids. The turbulent kinetic energy Ek is injected by stellar sources formed 
self-consistently in the simulation. In the ISM-like fluid, regimes of both forced and 
decaying turbulence coexist. In the active turbulent regions (forced regime), Ek is 
dissipated locally and efficiently. In the decaying regime (far from input sources), 
Ek(t) decays ~ (1 + t)~ - 8 . The residual turbulent motions may propagate dis- 
tances of the order of the observed disk height, suggesting that turbulence may be 
the responsible of vertical support and star formation self-regulation at the disk 
level, but not at the level of the whole cosmological halo, as would be required in 
some models of galaxy formation. 
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1. INTRODUCTION 

Modeling galaxy formation and evolution requires of a solid cosmological theoretical framework as much as 
an adequate model to describe the large-scale star formation (SF) cycle and its interplay with the interstellar 
medium (ISM). The dissipative properties of the ISM play a crucial role in this latter process. Globally, stellar 
radiation is mainly responsible for maintaining the temperature of the ISM in its various phases. However, 
thermal pressure is thought to be negligible for the global disk gas dynamics (e.g., Vazquez- Semadeni et al. 
1999; Cox, these proceedings). Nevertheless, stars are also sources of kinetic energy (Ek) deposition into the 
ISM, and as hydrodynamical simulations have shown, the dynamics of the gas in this case is deeply affected 
(e.g., Navarro & White 1993). Due to the large E k input and the high Reynolds number of the ISM plasma, 
turbulence is expected to develop, its pressure and dissipation being key ingredients in the "metabolism" of 
the disk stellar-gas system. 

Several disk galaxy evolution models (e.g., Firmani, Hernandez, & Gallagher 1996) are based on the idea 
that the intrinsic SF rate (SFR) is controlled by a balance within the vertical disk gas between the turbulent 
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energy input rate due to SF and the dissipation rate. The crucial parameter for the SFR and disk height is the 
turbulent dissipation timescale. 

The self-regulating SF mechanism has been also used in models of galaxy formation within the context of 
the hierarchical CDM-based scenario, but in this case it was applied to the large cosmological halo (White & 
Frenk 1991; Kauffmann, White, & Guiderdoni 1993; Cole et al. 1994; Somcrville & Primack 1999; van den 
Bosch 1999). In these models the feedback of the stars is assumed to efficiently reheat and drive back the disk 
gas into the dark matter halo, in such a way that the SFR efficiency is a strong function of the halo mass. Thus, 
a crucial question is whether the energy released by SNe and stars is able to not only maintain the warm and 
hot phases and the stirring of the ISM, but also to sustain a huge hot corona in quasi-hydrostatic equilibrium 
with the cosmological halo. 

It should be emphasized that the observed medium around the disks — diffuse ionized and high-velocity- 
dispcrsion HI gas, usually called the halo — is much more local than the hypothetical gas in virial equilibrium 
with the huge dark halo. We shall refer to the former as the the extraplanar medium, and to the latter as the 
intrahalo medium. It is still not at all clear how to explain the observed extraplanar medium, particularly the 
ionized gas (see a recent review by Mac Low 1999). This calls into question the possibility that ionizing sources 
from the disk (mainly massive OB stars) are able to sustain the extended intrahalo medium. A possibility is 
that this gas is heated by turbulence from the disk. However, this question again depends on the ability of 
the turbulent ISM to dissipate its Ek- Avila-Reese & Vazquez- Semadeni (in preparation; hereafter AV) have 
studied the dissipative properties of compressible MHD fluids that resemble the ISM. Here we briefly report 
their main results and remark the implications on the aforementioned questions. 

2. THE METHOD 

AV have used numerical 2D MHD simulations of self-gravitating turbulent compressible fluids that include 
terms for radiative cooling, heating, rotation and stellar energy injection (Vazqucz-Semadeni, Passot, & Pouquct 
1995,1996; Passot, Vazquez-Semadeni, & Pouquet 1996). The parameters were chosen in such a way the 
simulations resemble the ISM in the plane of the Galaxy at the 1 kpc scale. 

Previous simulations on dissipation in compressible MHD fluids were focused to study molecular clouds 
(Mac Low et al. 1998; Stone, Ostriker, & Gammie 1998; Padoan & Nordlund 1999; Mac Low 1999). In those 
works where the forced case was studied, the turbulence was driven in Fourier space by large-scale random 
velocity perturbations whose amplitudes are selected as to maintain Ek constant in time. As a result, Ek 
is injected everywhere in space (a "ubiquitous" injection). Instead, in the ISM the stellar input sources are 
pointlike and their spheres of direct influence are comparatively small w.r.t. to typical scales of the global 
ISM. In the simulations of AV, an "energy input source" is turned on at grid point x whenever p(x) > p c , and 
V • u{x) < 0. Once SF has turned on at a given grid point, it stays on for a time interval At s , during which the 
gas receives an acceleration a directed radially away from this point. The input sources are spatially extended 
by convolving their spatial distribution with a Gaussian of width A/. For the turbulent fluid, A/ is the forcing 
scale. At this scale the acceleration a produces a velocity difference around the "star" vj ~ 2aAt s (the velocity 
at which turbulence is forced at the A/ scale). Both A/ and vj are free parameters. 

3. DISSIPATION IN DRIVEN AND DECAYING REGIMES 

For ISM simulations (128 2 ) with driven turbulence, AV find that the behaviour with time of the Ek dissi- 
pation rate, E%, is similar to that of the Ek injection rate, Ek- This means that Ek is dissipated locally, near 
the input sources. For various simulations varying A/ and Vf, it was found that the the dissipation timescale 
is given by t d w 1.5 - 3.0 10 7 ( ^/gp^-i ), i- e -> t d is proportional to Xf/v f . 

Due to the locality and discretness of the energy input sources, most of the volume actually is occupied by 
a turbulent flow in a decaying regime. Thus, one may say that in the same fluid "active" turbulent regions, 
where the turbulence driven by small non-ubiquitous input sources is locally dissipated, coexist with extended 
regions with a "residual" turbulence in decaying regime and characterized by w rms , where u rms <C Vf. In order 
to study the decaying regime, SF was turned off in the simulations after some time t 3> td- It was found that 
Ek decays as (1 + t)~ n with n ~ 0.8, in good agreement with previous studies for isothermal fluids (Mac Low et 
al. 1998, Stone et al. 1998). A typical decaying timescale, idec, may be defined as the time at which the initial 



TURBULENT DISSIPATION IN THE ISM 



3 



Ek has decreased by a factor 2. From our simulations, tdcc « 1.7 x 10 7 years, which is in agreement with 
in the driven turbulence. With these timescales, "residual" turbulent motions propagating at roughly 10 km/s 
would attain typical distances of approximately 200 pc. It was also suggested, on dimmensional arguments, 
that Ek and v Ims will decay with distance £ as £~ 2m and £~ m , respectively, with m = n/(2 — n). 



4. CONCLUSIONS AND IMPLICATIONS 



• Localized, discrete forcing at small scales gives rise to the coexistence of both forced and decaying turbu- 
lence regimes in the same flow (ISM). 

• The turbulent Ek near the "active" turbulent regions is dissipated locally and efficiently. The global 
dissipation timescale td is proportional to A/ jvf. For reasonable values of A/ and Vf (which produce u rms ~ 10 
km/s), td is of the order of a few 10 7 years. Far from the sources, the "residual" ISM turbulence decays as 
Ek(t) oc (1 + t)~ - 8 . The characteristic decay time is again a few 10 7 years, and for v rms ~ 10 km/s, the 
turbulent motions reach distances of ~ 200 pc. 

• Turbulent motions produced in the disk plane will propagate up to distances of the order of the gaseous 
disk height. Therefore, models of galaxy evolution where this height is determined by an energy balance that 
self-regulates SF in the ISM appear viable. However, our results pose a serious difficulty for models of galaxy 
formation where the turbulent Ek injected by SNe is thought to be able to reheat and drive back the gas from 
the disk into the intrahalo medium in such a way that the SF is self-regulated at the level of the cosmological 
halo. Nevertheless, for non-stationary runaway SF (starbursts), most of the superbubbles might be able to 
blowout of the disk, as required for expelling large amounts of gas and energy into the dark matter halo. 
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